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A low-energy effective Yang-Mills theory for quark and gluon confinement 
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We derive a gauge-invariant low-energy effective model of the Yang-Mills theory. We find that the 
effective gluon propagator belongs to the Gribov-Stingl type and agrees with it when a mass term 
which breaks nilpotency of the BRST symmetry is included. We show that the effective model with 
gluon propagator of the Gribov-Stingl type exhibits both quark and gluon confinement: the Wilson 
loop average has the area law and the Schwinger function violates reflection positivity. However, we 
argue that both quark and gluon confinement can be obtained even in the absence of such a mass 
term. 
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It is well-known that the area law of the Wilson loop 
average is a gauge-invariant criterion for quark confine- 
ment. However, a gauge- invariant criterion for gluon con- 
finement and color confinement is not yet achieved. In 
recent several years, nevertheless, great endeavors have 
been made to clarify the deep infrared behavior of prop- 
agators for gluon and the Faddeev-Popov (FP) ghost in 
specific gauges, e.g., Landau, Coulomb and Maximally 
Abelian (MA) gauges. This research is motivated from a 
hope that color confinement might be attributed to the 
deep infrared behavior of the gluon and ghost propaga- 
tors 

In the most common Landau gauge, especially, it is 
still under debate to discriminate two different types of 
propagators, i.e., scaling 0] (an infrared suppressed gluon 
propagator with a finite or even vanishing dressing func- 
tion at zero momentum, and a ghost propagator more 
divergent in the infrared than its tree-level counterpart) 
and decoupling Q (an infrared finite gluon propagator 
and an ghost propagator with a finite dressing function 
at zero momentum) The scaling solution is in accordance 
with the Kugo-Ojima (KO) color confinement criterion 
1] and the Gribov-Zwanziger (GZ) confinement scenario 
2|. On the other hand, the decoupling solution is sup- 
ported by recent results of numerical simulations on the 
lattice with very large volumes See e.g. [1,0] on the 
present status of development. 

It is demonstrated [8[ that there is a one-parameter 
family of solutions for the ghost and gluon propagators 
of Landau gauge Yang-Mills theory and that it is only 
a matter of infrared boundary conditions whether in- 
frared scaling or decoupling occurs. Here the scaling so- 
lution is the only one member of this family that satisfies 
the KO/GZ property with a globally well-defined BRST 
charge. The remaining solutions are of a decoupling type 
and cannot maintain global color symmetry and Becchi- 
Rouet-Stora-Tyutin (BRST) symmetry simultaneously. 
Moreover, both type of solutions violates the reflection 
positivity which is a necessary condition for gluon con- 
finement, suggesting that neither type of solution can 
be associated with a massive gluon characterized by a 
gauge-independent pole mass. 

On the other hand, both type of solutions in the Lan- 
dau gauge has been shown to satisfy quark confinement 



criterion However, this result is valid for non-zero 
temperature T below the deconfinement temperature T c 
(0 < T < T c ), since vanishing Polyakov loop average at 
finite temperature was used as a gauge-invariant criterion 
for quark confinement [l(J EH • 

In this paper we examine whether there is a specific in- 
frared behavior of the gluon propagator which is compat- 
ible with both quark confinement and gluon confinement 
at zero temperature, irrespective of the gauge choice. For 
this purpose, we derive a gauge-invariant low-energy ef- 
fective model of the Yang-Mills theory at zero temper- 
ature in the gauge-independent manner. We show that 
the resulting low-energy effective gluon propagator be- 
longs to the Gribov-Stingl type [12| in the low-energy 
region. In the MA gauge, especially, the effective model 
is confining in the sense that the Wilson loop average 
has the area law and that the gluon Schwinger function 
violates the reflection positivity. In our model, an effec- 
tive gluon propagator agrees with the Gribov-Stingl form 
only when one includes a certain mass term violating the 
nilpotent BRST symmetry. However, we argue that such 
a mass term is not necessarily indispensable to obtain 
quark and gluon confinement simultaneously, since both 
the area law and positivity violation can be obtained even 
in the absence of such a mass term. 

This paper is organized as follows. 

(Step 1) [Reformulating the Yang-Mills theory in terms 
of new variables] In a path-integral quantization for the 
Yang-Mills theory, we decompose the Yang-Mills field 
s&n[x) into two pieces ^(x) and ^(x), i.e., stf^ix) — 
~ffi(x) + SEpix), and rewrite the action Sym[^] and the 
integration measure [d&f\ in terms of new variables re- 
lated to ~y^{x) and Sf^lx), according to filj - fi1| and fiH - 

(Step 2) [Deriving an effective model by eliminating 
high-energy modes] We integrate out field as the 
high-energy mode (p 2 > M 2 ) with a certain mass scale M 
of the field si^. Therefore, the resulting model S'ymM 
is written in terms of Y^{x), and is identified with a 
low-energy effective model for describing the low-energy 
regime p 2 < M 2 . A physical reasoning behind this step is 
explained below. The full gauge invariance of the original 
Yang-Mills theory SVm is retained also for SyM IX] ■ 

However, from the physical viewpoint of clarifying 
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what is the mechanism for confinement, we modify the 
step 1 and step 2 as follows. 

(Step 1', 2') We introduce an antisymmetric tensor 
field *B 11U of rank 2 (20l - [22l |. which is interpreted as a 
composite field of the Yang-Mills field. Then we repeat 
the same procedures as before to obtain the effective 
model SymIX, B] by integrating out 2£ field. 

By integrating out the B field in Sy^Y, B], an- 
other effective theory Sym[^] is obtained in the gauge- 
independent manner. 

(Step 3) [Converting the Wilson loop to the surface- 
integral] In the new formulation using new variables, we 
can exactly rewrite the Wilson loop operator Wc[^] 
originally defined in terms of srf^ix) by making use of 
alone without any reference to ^(x), according 
to dl^. This fact suggests that S^f] is suitable as 
a low-energy effective model for quark confinement. 

Up to now, all the results are gauge independent. In 
what follows, we choose a gauge to simplify the calcula- 
tions. 

(Step 4) [Calculating the Wilson loop average to show 
area law: quark confinement] The Wilson loop average 
(W / c[^'])ym, i-e., vacuum expectation value of the Wil- 
son loop operator Wc{^\ is evaluated by using the ef- 
fective model Sf M [f] as (WcKDym- We show that the 
Wilson loop average has the area law for sufficiently large 
loop C, leading to the non- vanishing string tension a in 
the linear part oR for the static quark- antiquark poten- 
tial V(R). 

Although the area law of the Wilson loop average is 
obtained also in the original model S Y m \V\ , the modified 
model Sym[Y,B] h &s the advantages: 

(Step 5) [Calculating the Schwinger function to show 
positivity violation: gluon confinement] The effective 
gluon propagator for "f of another effective theory 
SymI^I (obtained by integrating out the B field in 
Sf M [Y,B}) has the Gribov-Stingl form [13 • The 
Schwinger function calculated from the effective gluon 
propagator of the Gribov-Stingl form [l2j exhibits posi- 
tivity violation suggesting gluon confinement [l|, [|[ . 

Thus the derived effective model exhibits both quark 
confinement (area law) and gluon confinement (positivity 
violation). 

In this paper, we consider only the SU(2) gauge group 
[TUlH and the extension to SU(N) based on pj| will be 
given in a subsequent paper. 

(Step I) The explicit transformation from the original 
variables sd^ to the new variables V^, SC^ are given by 

Vn(x) = c^,{x)n{x) + ig~ 1 [n(x),d IJ ,n(x)], 

c^x) := ^(x) ■ n(x), 
X^x) = ig- l [D^n(x),n{x)]. (1) 

Here n(x) is the Lie-algebra sw(2)-valued field n(x) — 
n A (x)T A (A = 1,2,3) with a unit length, i.e., 
n A (x)n A (x) = 1. The so-called color direction field n 
must be obtained in advance as a functional of the orig- 



inal variable ^by solving the reduction condition [l6j . 
e.g., [n(x),D^]D tl [£/]n(x)] = 0. 

The new variable ~f)j,(x) as an su(2)-valued field 
f^x) = 'f A {x)T A (A = I, 2, 3) is constructed so that 

(i) has the same gauge transformation as the 
original field i.e., V^x) -> n(x)i / f ,(x)n(xy + 
ig~ 1 £l(x)d^Vl(xy and hence its field strength ^^[f] := 
d»% - - ig[y»,%] transforms as ^[Y^x) -> 
n(x)^ u [fMxy, and 

(ii) &p V \y\ is proportional to n, i.e., ,^^ v \f\{x) := 
n{x)G iiV {x). 

Consequently, G^ v — n ■ & tiv \¥\ is gauge-invariant, 
since the field n is constructed so that it transforms 
as n(x) — » D.(x)n(x)D,(xy . Remarkably, G^ LV has the 
same form as the 't Hooft-Polyakov tensor for magnetic 
monopole: 



G ^Lii — d^Cii 



ig 1 n- [d^dufi]. 



(2) 



(Step 1') We can introduce a gauge-invariant antisym- 
metric tensor field (* B) U y of rank 2 by inserting a unity 
into the path-integral |20l - t22j |: 



d D xl{(*B)^ 



1=1 VB exp 



- {an ■ ^[f] - pn ■ ig[^, ^])} 2 J , (3) 

where * is the Hodge dual operation. Here (too many) 
parameters 7, a, /3 are introduced to see effects of each 
term. When (3 = 7 _1 = G and a = 0, indeed, (*-B) Ml / 
is regarded as a collective field for the composite opera- 
tor n ■ ig\S£ il , with the propagator G obtainable in a 
self-consistent way [26| according to the Wilsonian renor- 
malization group (RG) |27j |. 

Then the Euclidean Yang-Mills Lagrangian is rewritten 
and modified into 



^ym\Y, 3£ , B] 



G\v + —{*B)^ V '—(* B)^Gf Jiiy 



I + 7a 2 

1 „„..„ „„ D 1+7/3 2 ,. 



2" 



(ig{^^})\ (4) 



where we have defined 

^S AB 8, V + (2 + ia p)ge ABC n c G^ 



S 



AB 



1 fige A ^n^{*B)^ 
(D p [-f]D p [Y]) AB , 



(5) 



with the covariant derivative in the adjoint represen- 
tation with ^ := yfTc, (T C ) AB = if ACB : D AB := 
d»6 AB - gf ABC YC = [d M l - ig^] AB . 

[On the effect and the role of the gluon mass term] The 
gluon "mass term" for the S£ field, 



(6) 
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is gauge invariant in the new formulation [181 ]. Therefore, 
we can include this mass term in calculating the low- 
energy effective action. But we do not introduce this 
mass term explicitly. On the other hand, the inclusion of 
the gluon mass term for the ~f field, 



1 2v/2 1 2 2 

2 ii 2 



(7) 



breaks gauge invariance and BRST invariance after tak- 
ing specific gauges. However, we can modify the BRST 
such that the modified BRST is a symmetry of the Yang- 
Mills theory with the mass term at the cost of nilpotency. 

(Step2') We identify 5£p with the "high-energy" mode 
in the range p 2 € [M 2 , A 2 ] and proceed to integrate out 
the "high-energy" modes X^. Here M is the infrared 
(IR) cutoff and A is the ultraviolet (UV) cutoff as the 
initial value for the Wilsonian RG. In the derivation of 
our effective model, we neglect quartic self-interactions 
among S£ p , i.e., (ig[S^ p , ^]) 2 ■ 1 In these approxima- 
tions, we can integrate out SC^ by the Gaussian integra- 
tion and obtain a gauge-invariant low-energy effective 
action SymI^i B] without mass terms <J6]) , (JTJ) : 



rl 



W r 2 
1 ^P* 



■m%,- 



'(* B) pa G pa 



-lndetQ^f -In det S" 45 , 



(9) 



where J = J d 4 x, the functional logarithmic determinant 
\ In det Q AB comes from integrating out the S€ field, and 
the last term comes from the FP-like determinant term 
17] associated with the reduction condition 16]: 



In det Q£f - In det S AB 



g 2 \n± 



M 2 



(4tt) 2 
2 1 1 



1 G 2 



--{{2 + ia p)G pa - 1 p{*B) pa Y 
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(4tt) 2 M 2 6 
0(d 4 /M 4 ). 



(d x {{2 + 1 ap)G pa ~ 1 P(*B) pa }) 



(10) 

The gauge fixing is unnecessary in this calculation. In- 
deed, the resulting effective action © with (JTUJ) is man- 
ifestly gauge invariant. This is one of main results. 



1 However, we can take into account an effect coming from the 
quartic interaction, which influences our effective model. In 
fact, it is shown 18, 33] that the quartic gluon interaction 
among &ffi gluons can induce a contribution to 
the mass term I© \M 2 X 2 through a vacuum condensation of 
"mass dimension-2" (the BRST-invariant version was proposed 
in 

(%?{x)X?(x)) (8) 

which leads to the mass term I© with M 2 ~ (%®{x)%®(x)) 
up to a numerical factor. This result is easily understood by a 
Hartree-Fock argument. This effect is included in the heat kernel 
calculation through the infrared regularization. 



The correct RG /3-function at the one-loop level 13(g) := 
M M££) = -b l9 3 + 0(g 5 ), bi = f /(4tt) 2 is reproduced in 
a gauge invariant way when ja(3 = which follows from 
e.g. a = (mentioned above) or 7 = (in the case of no 
B„ v field). 

To obtain (TTU1) . we used the heat kernel to calculate 
the regularized logarithmic determinant. Instead of using 
the standard regulator function Rm of the functional RG 
approach [27| . we restrict the integration range of r to 
t G [1/A 2 , 1/M 2 ], which corresponds to the momentum- 
shell integration p 2 g [M 2 ,A 2 ] 



In det = - 



x lim — 

s->o ds 



,2s 



r( s ) 



1/M 2 



dTT S 



1/A 2 



x tr(a;|e~ TO |a;) 



(11) 



where tr denotes the trace over Lorentz indices and group 
indices and /x is the renormalization scale. The limit A — > 
00 should be understood in what follows. These results 
should be compared with previous works [13, HH, [H, US] ■ 
We can show that the mass term © plays the same role 
as the IR regulator mentioned above, see [22[ . 

(Step 3) We use a non-Abelian Stokes theorem [23M25| 
to rewrite a non-Abelian Wilson loop operator 



WcH :=tr 



exp 



/g (f> dx^stf^x) 



(12) 



into the area-integral over the surface £ (9E = C): 

1 



WcW] 



d/i S (£)exp 



G 



(13) 



where the product measure dfiY.{C) '■= ELes ^fe) ^ s 
defined with an invariant measure dfi on SU (2) normal- 
ized as J d/j,(£ x ) = 1, € SU(2). In the two-form G := 
\G llu (x)dx^ A dx v , Gpy agrees with the field strength 
@ under the identification of the color field n(x) with 
a normalized traceless field n(x) := 6e(0"3/2)£!- See also 

(Step 4) We proceed to evaluate the Wilson loop av- 
erage W(C) = (Wc[^])ym by using the effective action 
Sf M [Y, B], i.e., (WcM)ym ^ {WcW\)f M with the aid 

of CESD. 

We demonstrate that the simplest way to obtain 
the area law is to use the low-energy effective action 
S'ym I c > B] retained up to terms quadratic and bilinear 
in c and B : 



S- 



off 
YM 



z,B] 



- 7a 



G: 



+ 



+ 



g 2 ln± 



A/ 2 



1, 



"2~(* B) pa Gpa 



(4tt) 2 

_£ i_I 

(4tt) 2 M 2 6 
1 



ia p)G pa - iP{*B) pa ] 2 
{d x [{2 + 1 ap)G pa - 1 p(;B) pa \j 2 



+ / * m < + o(a 4 /M 4 ), 



(14) 
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if we include an optional mass term (J7J . 2 

To obtain the propagator or correlation functions, we 
need to fix the gauge. For instance, in the Landau 
gauge, d^s^n = 0, correlation functions for new variables 
have been computed on a lattice by numerical simula- 
tions using the Monte-Carlo method in [36] based on 
pjil HH. This justifies the identification of <^ as the 
high-energy mode negligible in the low-energy regime be- 
low M ~ 1.2GeV. 

In what follows, we take the unitary-like gauge 



n A (x) = S A 3, 



(15) 



which reproduces the same effect as taking the MA 
gauge [37] in the original Yang-Mills theory. In this 
gauge, %Mx) reduces to the off-diagonal component 
A^(x) (a = 1,2), while ~f A {x) reduces to the diag- 
onal one A^(x) = a^x), i.e., % A {x) = ^(x)5 Aa , 
Y l f(x) = £/^(x)Sa3 — c^(x)5a3- In this gauge, the field 
strength reads 

G^(x) -> F tiv {x) := d fl c v (x)-d u c fl (x), F = dc, (16) 

where d denotes the exterior differential. The gauge (fT"5)) 
forces the color field at each spacetime point to take the 
same direction by gauge rotations. Hence the field c 
contains singularities (of hedge-hog type) similar to the 
Dirac magnetic monopole after taking the gauge (fT5|) . 
Therefore, dF = ddc ^ 0, even if F = dc. If we do 
not fix the gauge, such a contribution is contained also 
in the part ig~ 1 n ■ [d^n, d v n] to make a g auge-invariant 
combination given by © , see [H], Ha] . 

By integrating out the B field, we obtain the effective 
action SyL [c] . Then we find that the effective propagator 



has the Gribov-Stingl form: 



1 + d lP 2 



c + cip 2 + C2P 



(17) 



TO 2 , Cl 



•< „ g 1 iyo 1 



where Co 

jaf3) 2 + (1 + 7a 2 h^ 2 + 2(2 + ■ya/3)ja/3]/3, and di = 
~*f~ r£pp jj2- The precise values of the parameters 
m, 7iO, P and M are to be determined by the functional 
RG [13] following [ll|, which is a subject of future study. 



2 The optional mass term is not contained in the original ac- 
tion of Yang-Mills theory. But, there is a possibility that such a 
mass term could be induced in a non-perturbative manner in 
the effective theory for the deep infrared region. In the full 
non-perturbative treatment, it is important to avoid the Gri- 
bov copies to achieve complete gauge fixing. We observe that, in 
the modified version of the GZ model Q| called the refined GZ 
[3 ll | . the similar effect to our optional mass term is generated by 
the horizon term which plays the role of restricting the range of 
the functional integral to the first Gribov region. The horizon 
term breaks the nilpotent BRST symmetry, so does the optional 
mass term 10- This issue is skipped here. 



In the unitary-like gauge (fTS"]) the Wilson loop operator 
is reduced to 



W C [F] = exp 









= exp 


. 2 Js : as=c _ 





*s^e s ,F) 



where Ox; is the vorticity tensor defined by Q^(x) — 
J s d 2 S^ v {x(a))5 D (x— ir(<7)), which has the support on the 
surface £ whose boundary is the loop C. Here (•, •) is the 
L 2 inner product for two differential forms: (0-£,F) — 
J d D x\Q l ^{x)F tlv {x). By integrating out B, we obtain 

the effective model SymM = 5 ( c ' ^<x^) 



FF 



in terms of c or F. Then the Wilson loop average W(C) 
is evaluated by integrating out F = dc: 



W(C) = exp 



1 



3 2 (e, 



?FF 



(19) 



where &pp — A3) cc and its Fourier transform &pp{p) = 
P 2 &cc(p)- For concreteness, we choose Os for a pla- 
nar surface bounded by a rectangular loop C with side 
lengths T and R in the X3 — xa plane. Then we find that 
the Wilson loop average has the area law for large R 



W{C) ~ exp[-CTi?T], 
with the string tension given by the formula: 



(20) 



p{+Pi<M 2 



dpidp 2 % 

(27T) 2 



?FF 



(pi,P2,0,0)>0, (21) 



where the momentum integration is restricted to the 
two-dimensional space (the dimensional reduction by two 
[Hj|) and is cutoff at the upper limit M. A positive and fi- 
nite string tension < a < 00 follows from the condition 
of no real poles in the effective gluon propagator 2# cc (p) 
in the Euclidean region, < 3>pp(p) — p 2( 2> cc {p) < 00, 
which is connected to gluon confinement below. This is 
another of main results. 

According to numerical simulations in MA gauge [39r - 
l4lj . the diagonal gluo n propagator is well fitted to the 
form (HU): e.g. [41] give c = 0.064(2) GeV 2 , c x = 
0.125(9), c 2 = 0.197(9)GeV~ 2 , d x = 0.13(l)GeV -2 , and 
M ~ 0.97GeV, where M is the mass of off-diagonal glu- 
ons obtained in the MA gauge. This value of M is a little 
bit smaller than the values of other groups [39l . |40| . This 
indeed leads to a good estimate for the string tension a ~ 
(0.5GeV) 2 according to for a(/i) = 3 2 (^)/(4tt) ~ 1.0 
at /j, = M. The next task is to study how the results are 
sensitive to the deep infrared behavior of the diagonal 
gluon propagator (|17l) and the actual value of M for the 
off-diagonal gluon propagator. 

The Gribov-Stingl form is obtained only when Co 7^ 
(i.e., to 7^ 0) and d\ 7^ (B^ is included). Even in the 
limit to 2 — > (co — > 0), the area law survives according 
to ([2Tj) , provided that @pp(p) remains positive and finite: 
$>pp(p) — > C ^ + C2P I ' wnne ^cc(p) behaves unexpectedly as 

@cc(p) -> pil^+czp*) - Hence, we argue that it does not 
matter to quark confinement whether m — or to 7^ 0. 
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(Step5) The positivity violation is examined. In the 
case of c 2 = 0, there is no positivity violation, as far as 
co/cx > 0. In the case of c 2 ^ 0, @ cc {p) has a pair of com- 
plex conjugate poles at p 2 — z and p 2 — z* , z := x + iy, 

x := -ci/(2c 2 ), y := \fc /c 2 - (ci/(2c 2 )) 2 . We find that 

the Schwinger function A(t) := f^e ip4t # cc (p = 

0,^4) is oscillatory in i and is negative over finite in- 
tervals in the Euclidean time t > 0: 

aw = 2c2N 3 /2 1 sin(2y) e^' 1/2 -Ht N v^ os ,-, ) 

+ di|z|cos(i|z| 1/2 cos<^ + < / 3)], (22) 

where z = \z\e 2lip with |z| = (c /c 2 ) 1 ^ 2 , cos(2(^) = 
-V / c 2 /(4c c 2 ), and sin(2<^) = y/l - c 2 /(4c c 2 ). There- 
fore, the reflection positivity is violated for the gluon 

c 2 

propagator p7|) . as long as < 4 ^ < 1, irrespective of 
di. When Co = (or m = 0), 

A(*) = -^ * /^(i-^VVI. (23) 
2ci ZCx V c l V c 2 / 



Hence, the special case Co = also violates the positivity, 
if Ci > and c 2 > 0. Thus the diagonal gluon in the MA 
gauge can be confined. 

In summary, we have derived a novel low-energy ef- 
fective model of the SU(2) Yang-Mills theory without 
fixing the original gauge symmetry. It is remarkable that 
the effective model respects the SU(2) gauge invariance 
of the original Yang-Mills theory, which allows one to 
take any gauge fixing in computing physical quantities 
of interest in the low-energy region. The effective gluon 
propagator belongs to the Gribov-Stingl form. In MA 
gauge, the model exhibits both quark confinement and 
gluon confinement simultaneously in the sense that the 
Wilson loop average satisfies the area law (i.e., the linear 
quark-antiquark potential) and that the Schwinger func- 
tion violates reflection positivity. More results and full 
details will be given in a subsequent paper. 
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